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We propose a new system suitable for studying analogue gravity effects, consisting of a gas ffow¬ 
ing in a duct with a compliant wall. Effective transonic flows are obtained from uniform, low 
Mach number flows through the reduction of the one-dimensional speed of sound induced by the 
wall compliance. We show that the modified equation for linear perturbations can be written in a 
Hamiltonian form. We perform a one-dimensional reduction consistent with the canonical formu¬ 
lation, and deduce the analogue metric along with the first dispersive term. In a weak dispersive 
regime, the spectrum emitted from a sonic horizon is numerically shown to be Planckian, and with 
a temperature fixed by the analogue surface gravity. 

PACS numbers: 04.62.+v, 04.70.Dy, 43.20.+g, 43.20.Wd 


Engineering flows that are transonic and regular offers 
the possibility to test well known predictions concern¬ 
ing astrophysical black holes [1]. Of particular interest 
is Hawking’s discovery that black holes should sponta¬ 
neously emit a steady thermal flux [2]. Although this 
effect was originally phrased in the context of quantum 
relativistic fields, it rests on the anomalous mode mix¬ 
ing occurring near the black hole horizon [3]. This mix¬ 
ing, which is stationary and conserves the wave energy, 
is said anomalous as it leads to a mode amplification and 
involves negative energy waves. Because of the precise 
analogy between the equation governing sound propaga¬ 
tion and that used by Hawking, these key elements are 
recovered in a stationary transonic flow. Indeed, in the 
acoustic approximation, for long wavelengths, the mode 
mixing possesses the main properties of the one respon¬ 
sible for the Hawking effect [4, 5]. 

To complete the comparison, one should take into ac¬ 
count the dispersive properties of sound waves, which 
have no counterpart in general relativity. (Note however 
that dispersive terms appear in certain theories of mod¬ 
ified gravity where Lorentz invariance is broken [6-8].) 
Analytical and numerical studies have established that 
the correspondence is quantitatively preserved provided 
the two relevant scales are well separated [4, 5, 9-12], 
namely when the dispersive length is sufficiently smaller 
than the typical length scale associated with the inhomo¬ 
geneity of the flow (which then plays the role of the in¬ 
verse surface gravity of the black hole). Therefore, there 
is no conceptual obstacle preventing to test the Hawk¬ 
ing prediction by observing the mode mixing across a 
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EIG. 1. Schematic drawing of the configuration. At y = 0, 
the wall is rigid. At y = 1, the compliant wall is made of a 
succession of tiny tubes of height b{x). As explained in the 
text, the effective sound speed is a decreasing function of b. 
In the present prohle b(x) and for a given uniform flow with 
Mach number M < 1, an effective supersonic region can be 
created in the right region. 

sonic horizon. In practice, the difficulty is to find appro¬ 
priate set-ups. Many have been proposed, involving for 
instance ultra cold atomic clouds [13], surface waves in 
flumes [14], and light in non-linear media [15]. Recently, 
the first experiments have been carried out [16-19]. 

In this Letter, we propose a new framework which is 
a variant of the original one [1, 4]. Its main interest re¬ 
sides in obtaining a large reduction of the low frequency 
one-dimensional sound speed in a duct so that a station¬ 
ary flow with a uniform low Mach number M possesses 
a horizon. In realistic settings, M could be close to 0.3. 
The reduction of the effective sound speed is achieved 
by means of a compliant wall composed of thin tubes 
which modify the upper boundary condition, see Fig. 1. 
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Interestingly, for near-critical flows, the wave equation 
follows from a well-defined action principle. From this, 
we derive the acoustic black hole metric in the long wave¬ 
length regime, and a conserved norm which establishes 
the anomalous character of the mode mixing. At the end 
of this Letter we briefly present the practical advantages 
of this proposal. 

The model - We consider the propagation of sound 
waves in a 2-dimensional horizontal channel of uniform 
height H. x* denotes the cartesian horizontal coordinate, 
the vertical one, and t* the time. We assume the flow 
of air is uniform, with a horizontal velocity Uq. Denoting 
by Co the sound speed, po the air density, v* the velocity 
perturbation, and p* the pressure perturbation, the time 
evolution is given by 

= -poV* ■ V*, poDtv* = -W, (1) 

where = dt* + UqOx* is the convective derivative. We 
define dimensionless quantities as x = x*/iL, p = p*/iL, 
t = t*co/H, V = v*/co, p = p*/{poCo), and M = Uq/cq. 
The potential (j) gives the velocity by v = V0, and the 
pressure hy p = —Dtcj). It obeys 

D^^-{dl + dl)4> = o. ( 2 ) 

At the lower wall, the impenetrability condition is sim¬ 
ply dy(j){t^x^y = 0) = 0, see Fig. 1. At the upper wall 
p = 1, the continuity of the displacement and pressure 
gives rise to a non local expression in time, see [20-22] 
for details. However, for near-critical flows and small 
frequencies, it can be written as 

dycj) + Dt{h{x)Dt(t)) = 0 at p = 1, (3) 


We then adjust the flow speed M to obtain 
1 1 


cs{h) = 




< \M\ < 


b2 


= cs{b2). ( 8 ) 


Since the background flow is stationary, we shall work 
with (complex) stationary waves: 


f 

= / 


df . 

(9) 


Wave equation - Since the height H of the duct is 
much smaller than typical longitudinal wavelengths, we 
expect that stationary waves obey an effective 1 dimen¬ 
sional equation in x, as it is the case in elongated atomic 
Bose condensates [23], and in flumes [14, 24]. To obtain 
such a reduction is non trivial as the x-dependence of b{x) 
prevents us from factorizing out a p-dependent factor. To 
proceed, and to make contact with the above Refs., it is 
useful to exploit the fact that Eq. (2) and Eq. (3) can be 
derived from the following action 

S = tJdtJdxJ dyC 

C = {Dt<pf - (V</.)2 + S{y - l)b{x){Dt4>f. (10) 


We notice that the non-trivial condition of Eq. (3) is in¬ 
corporated by the above boundary term. Introducing the 
conjugate momentum 7 r(t, x, p) = dC/d{dt(t)) one obtains 
the Hamiltonian H by the usual Legendre transformation. 
In addition, as in [4, 24], the conserved inner product (-I*) 
is not positive definite, and has the Klein-Gordon form. 


which is second-order in dt. 

Eor a homogeneous stationary flow, we can look for 
solutions of the form pk oc cosh((a/cp)e*^^^“^^^\ Eq. (2) 
and the two boundary conditions respectively give 

a\ = and(a/ctanh((a/c) = bVt\^ (4) 

where = uj^ — Mk is the frequency in the co-moving 
frame. At low wave number, the dispersion relation reads 


nl = cUb)k^-kyAl + 0{k% (5) 


where 


cl{b) = 


- 


3(1 + 6 )^ 
62 


(6) 


One sees the important effect of the boundary condition 
of Eq. (3): the low frequency group velocity with respect 
to the fluid (= 9/^0) is reduced by the compliant wall. 
One also sees that the dispersive length I/A 5 given by 
the quartic term vanishes in the limit 6 ^ 0 . 

To obtain flows crossing the effective sound speed, we 
make b vary with x, see Eig. 1 . We call bi > 62 its 
asymptotic values, and its typical variation length. 
We choose the following form for b{x): 


b{x) 


bi -1- 62 


62 - bi 


2 


tanh 



( 7 ) 


/ -|-CO 1*1 

dx dy {Tr^(j )2 - (I>iTT 2 ) , (H) 

-00 J 0 

where 0i, 02 are two complex solutions of Eq. (2, 3), and 
7 ri, 7 r 2 their associated momenta. Eq. (11) is conserved 
in virtue of Hamilton’s equations. As for sound waves in 
other media, it identically vanishes for all real solutions. 
However it provides a key information when studying sta¬ 
tionary modes, namely the sign of their norm ( 0 a;, 0 a;). 
Indeed, as we shall see, for a fixed cj > 0, there will be 
both positive and negative norm modes. When the flow 
is stationary, for any complex solution 0 , the wave energy 
is conserved and related to Eq. (11) by 

H[2Re(0)] = (0,i9t0). (12) 

Moreover, when the flow is also asymptotically ho¬ 
mogeneous, for every asymptotic plane wave pk oc 
cosh{aky)e'^^^^~^^^\ the sign of H is that of uJu^k- This 
relation will allow us to identify the negative energy 
waves without ambiguity. 

We can now proceed following the hydrodynamic treat¬ 
ment of [24]. As a first step, it is useful to derive a (1+1)- 
dimensional equation from which an effective space-time 
metric can be read out. When the (adimensional) wave¬ 
length in the x direction is much larger than 1 , we can 
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assume that dy(j) is independent of As = 0 at 
^ = 0, we write the field as 

(j){x,y,t)^^{x,t)+y^'^{x,t). (13) 


Plugging this into the action Eq. (10) and varying it with 
respect to ^ and we get two coupled equations. Com¬ 
bining them, we obtain {O 2 + = 0, where On is a 

n^^-order operator in dt and The quadratic term is 
O 2 = d^F^^{x)dn, where 


F^^{x) = 


— '-s 


cl{h{x))-M^ M 


M 


(14) 


and Cg{b{x)) is given in Eq. (6). Up to a conformal factor, 
we obtain the d’Alembert equation in a two dimensional 
space-time with metric oc F^^. This metric has a 
Killing horizon where c|(x) — vanishes [25]. This 
correspondence with gravity relates the anomalous scat¬ 
tering described below to the Hawking effect. 

Contrary to what happens for sound waves in atomic 
BECs, or water waves in the incompressible limit, O 4 
also contains third and fourth derivatives in time. This 
prevents to apply the standard treatment on the sole field 
4>. However, the set of two coupled equations on (4>, T) 
is hamiltonian and can be used to study the scattering. 
Alternatively, one can work with the original model in 2-1- 
1 dimensions based on Eq. (10). We performed numerical 
simulations with both models and found similar results. 

Anomalous mode mixing - Since stationary waves with 
different frequencies uj do not mix, the scattering only 
concerns the discrete set of modes with the same uj. To 
characterize it, we identify its dimensionality and the 
norms of the various asymptotic modes for x Too. 


U) 



FIG. 2. Dispersion relation uj versus /c in a homogeneous 
subsonic flow. The blue, plain curves show the roots with 
positive comoving frequency H, and the red, dashed ones those 
with H < 0. The dotted black line represents uj = 0.4, and the 
dot-dashed one shows the value Ct;max of uj at which the two 
roots with /c < 0 merge. The parameters are 5 = 1, M = 0.4, 
and the effective sound velocity cs is equal to \l\f2 0.71. 

Eig. 2 shows the dispersion relation and the roots at 
fixed cj in a subsonic flow for M > 0. Eor this sign of M, 


the flow associated with Eq. (7) passes from supersonic 
to subsonic along the direction of the stream. It thus 
corresponds to a white-hole flow, as those studied in [17, 
24, 26]. ^ Eor definiteness, we discuss only the case uj > 
0. The same results are directly applicable to cj < 0 
after complex conjugation. In the subsonic region, on 
the right of the horizon, there exists a critical frequency 
<^max (close to 0.42 in the Eigure) at which two roots 
merge. Eor uj < cjmax, the dispersion relation has 4 real 
roots. ^ Eollowing [24], we call their wave-vectors 

5 and —k^^. The corresponding asymptotic 

modes are, respectively, and ((/:?UjJ)*. 

They are characterized by three important properties: 

1) in- or out-character: p]^ is incoming (it moves to¬ 
wards the horizon) while the three other modes are out¬ 
going (they move away from the horizon); 

2) energy sign: {p^^) carries a negative energy and 
a negative norm, see Eq. (12). (It has been complex 
conjugated so that p^^ is a positive norm mode). The 
three other modes have positive energy and norm; 

3) CO- or counter-propagating nature: (^co,out -g 
propagating (its group velocity in the frame of the fluid is 
positive) while the three others are counter-propagating. 
This separation is useful because only the latter are sig¬ 
nificantly mixed in a transonic flow [9, 10]. In effect, 
(^co,out essentially as a spectator. 

To get the S matrix, we need to identify the basis of in¬ 
coming (outgoing) modes which contain only 

one incoming (resp. outgoing) asymptotic plane wave. 
Eor UJ < cjniax, there are three modes, so the scatter¬ 
ing matrix has a size 3x3, and is an element of the 
Lie group U(2,1) since {p-ooT has a negative norm [9]. 
In this paper we focus on the incident mode for a 
white hole flow. It is a good candidate to probe the ana¬ 
logue Hawking effect, and was studied in hydrodynamic 
flows [16, 17, 26, 27]. Eor x ^ oc, it is a sum of four 
asymptotic modes 

lin ^ ,An , ^ ,^out , o , a ,^co,out 

9uj ^ F a^p^ + Puj • ( 13 ) 

In transonic flows, there is no transmitted wave [26]. 
thus vanishes for x -A- — oc. When working with asymp¬ 
totic modes of unit norm, the norm of plf evaluated at 
late time (in the sense of a broad wave packet) 

{N-A^ = \ap-\/3p + \AA, (16) 

must be exactly 1 because of the conservation of Eq. (11). 
(In stationary flows, (A^°^^)^ = 1 also expresses the con¬ 
servation of the wave energy, see Eq. (12), and that of 
the energy flux of Mohring [28].) The minus sign in 


^ For M < 0, one would describe a black hole flow. The forthcom¬ 
ing analysis applies by reversing the sign of velocities and the 
“in” or “out” character of the modes, see [9]. 

^ In the supersonic region, only two real wave-vectors remain: 
and They both describe incoming modes. 
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front of \Puj?‘ is the signature of an anomalous scatter¬ 
ing. It stems from the negative norm carried by (0°^) , 
see the above point 2. The coefficient (3^ thus mixes 
modes of opposite norms and energies. In quantum 
settings, \/3^\ would give the mean number of sponta¬ 
neously produced particles from amplifying vacuum fluc¬ 
tuations, that is, the Hawking radiation [2]. The gravi¬ 
tational analogy [1, 4] indicates that \/3ujf should follow 
a Planck law when dispersion effects (and grey body fac¬ 
tors [29, 30]) are negligible. Moreover, it predicts that the 
effective temperature T should be given by Th = /^/27r, 
where n is the surface gravity obtained from the ana¬ 
logue metric of Eq. (14). (We choose the units so that 
= 1. Th and n are thus both frequencies.) Using 
Eq. (6), one gets 


^ — ^x^s\cs=M — 2 ^x^\cs=M • 

Spectral analysis - We numerically solved the set of 
coupled equations on the (1 + l)-dimensional fields 4> and 
T, using the method of [26] adapted to the present case. 
The results concerning the incoming mode of Eq. (15) 
propagating in a transonic flow described by Eq. (7) are 
shown in Eig. 3. ^ We stopped the integration for uj 
slightly below the critical frequency c^max, where and 
both vanish. We tuned the various parameters (given 
in the caption of Eig. 3), so that the flow is near-critical: 
M/cs{bi) « 1.054 and M/cs{b 2 ) ~ 0.943. Using these 
parameters, one has cjmax ~ 0.0053, and k ~ 0.019 of the 
same order as the dispersive frequency scale A^cl evalu¬ 
ated at the horizon. This means that we worked just out¬ 
side the weak dispersive regime [11]. Yet, for frequencies 
up to Climax 5 follows rather well the Hawking predic¬ 

tion p = XjicAl^^ — 1), that is, a Planck law with Th 
given by /^/(27r), see Eq. (17). At low frequency, the rel¬ 
ative difference — 1 is of order 20% (when 

we used = 3, the difference reduced to about 0.3%, 
as expected since we were then in a weakly dispersive 
regime). Moreover, we see that the coefficient involv¬ 
ing the co-propagating mode can be safely neglected as 
\A^f remains smaller that 0.1%. In order to estimate the 
numerical errors, we also show the quantity (Y^^^)^ — 1 
which must be equal to zero, as explained below Eq. (16). 
In brief, the properties we obtain are in close agreement 
with those found in other media [9, 12, 23, 31]. 

Experimental aspects - The incident waves will be sent 
by a loudspeaker in a rectangular channel, and the scat¬ 
tered waves will be observed using two arrays of micro¬ 
phones, see Eig. 1. The frequency range will be chosen 
to fulfill the low frequency hypothesis used in Eq. (3). A 
stationary flow with a mean Mach number M ^ 0.3 will 
be provided in the channel. The compliant wall will be 
realized with a honeycomb structure. Its height H b{x) 
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FIG. 3. Plot of ^ \l3ujf (blue, plain) and \Auj\‘^ (orange, 
dashed) of Eq. (15) as functions of the frequency. The black, 
dot-dashed curve shows ^ for a Planck law at the 

Hawking temperature Th = k,/{ 27 t). The parameters are: 
M = 1/3, = 9,^2 = 7, and db = 1. The green, dotted line 

represents - 1 where is given in Eq. (16). Its 

non-vanishing value quantifies the numerical errors. 


will vary in such a way that the flow is transcritical, i.e., 
that M/cs crosses unity at some x, something which has 
not yet been reached in water tanks experiments [16, 17] 
aiming at detecting the analogue Hawking radiation. An¬ 
other important advantage is the possibility to send the 
three types of incident waves, and therefore to measure 
the 9 scattering coefficients, both in black hole and white 
hole flows. We suspect that working in these conditions 
can produce turbulence effects and whistling. To reduce 
these effects, the compliant wall will be covered by a wire 
gauze with a very low flow resistance. Special attention 
should also be devoted to dissipation. In spite of these 
experimental difficulties, the system seems to be a good 
candidate to probe the various aspects of the analogue 
Hawking radiation. 

Conclusions - We showed that a low Mach number uni¬ 
form flow of air in a tube with a compliant wall can pro¬ 
duce a sonic horizon by reducing the local effective one¬ 
dimensional speed of sound. Despite the unusual bound¬ 
ary condition at the compliant wall, the problem was 
phrased in a Hamiltonian formalism. Eor near-critical 
flows, a (1+1) dimensional-reduction was performed, ex¬ 
hibiting an analogue metric and the first dispersion terms 
while retaining the hamiltonian structure. The Hawk¬ 
ing spectrum was numerically recovered for sufficiently 
slowly varying profiles of the compliant wall. We thus 
hope to be able to verify that the norm of the anomalous 
coefficient grows as K,/{27ruj) for low frequency. 

Acknowledgements - We thank Scott Robertson for his 
remarks on a preliminary version of this letter. PE is 
grateful to the LPT for its hospitality during an intern¬ 
ship sponsored by the Ecole Normale superieure. 


3 When sending a localized wave-packet on a white hole horizon, observed at late times the formation of an undulation, see [22] . 
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Wave equation and small-frequeney approximation 

When considering a stationary mode oc the boundary condition at the compliant wall y = I reads [21] 


dy4>uj = - + Mdx) 


tan (ujb{x)) 


UJ 


{—iiO + Mdx) I 0a 


(SI) 


To avoid having to deal with time derivatives of high orders, we consider a small-frequency limit tan (u; 6 ) /cj ~ 6 , 
to obtain Eq. (3). It must be noted that this simplification does not affect the quadratic term in the dispersion 
relation, but changes the quartic term. The relative difference is of order (l — M/cs) when considering the counter- 
propagating modes. This approximation is thus valid for all the spanned range of cj G [0,co’max] provided M/cs{b{x)) 
remains everywhere close to unity. 

Wave equation and normalization in the (l-l-l)D model 
The two operators O 2 and O 4 appearing below Eq. (13) are given by 


© 2 = A {l + h{x))Dt-dl, 

A = ^ (a" - 52 ) (a (1 + 66(a;)) A - . (S2) 

Asymptotic modes are given by two-component plane waves: 

^x,t) = 

where T are defined in Eq. (13). The prefactors u and v are given by 

u 1(1 + 36) — k‘^ 

V ^ “3 (l + 6)fl2 ’ 

((1 + b)Mn + k) |m|^ + +(1 + 36) Mil + k) Re {uv*) + +(1 + 56) 

o o 

The first condition comes from the wave equation derived by varying the action Eq. (10) with respect to while the 
second one ensures that all the asymptotic modes have a unit norm (up to a minus sign). 

Wave-paeket and undulation 

In Eig. SI we show a space-time diagram of the perturbation obtained by sending a localized wave-packet on a white 
hole horizon. The important point is the appearance of a long-lasting undulation^ i.e., of a zero-frequency mode with 
a macroscopic amplitude. Its presence is due to the diverging character of ^ as l/uo for u; ^ 0 which 

amplifies the small-frequency components of the incident wave-packet. See Eq. (23) in [24] and Eq. (27) in [32] for 
studies of the same mechanism in related contexts. 


(53) 

(54) 

+ k) \v\^ = ± 1 . (S5) 
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FIG. SI. Space-time plot of the undulation obtained when sending a wave packet initially centered around x = 50 in a white 
hole flow. The sonic horizon is located around x = 0. As in the main text, the sign of M is positive and the subsonic region 
is on the right side of the horizon. The amplitude represented is that of F = c^dy = $ + |T. The top plot shows the 
time-evolution with a linear color scale, and the bottom one shows the late-time conhguration. 
















